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Abstract. We introduce a continuous data assimilation (downscaling) algo¬ 
rithm for the two-dimensional Navier-Stokes equations employing coarse mesh 
measurements of only one component of the velocity field. This algorithm can 
be implemented with a variety of finitely many observables: low Fourier modes, 
nodal values, finite volume averages, or finite elements. We provide conditions 
on the spatial resolution of the observed data, under the assumption that the 
observed data is free of noise, which are sufficient to show that the solution of 
the algorithm approaches, at an exponential rate asymptotically in time, to the 
unique exact unknown reference solution, of the 2D Navier-Stokes equations, 
associated with the observed (finite dimensional projection of) velocity. 
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1. Introduction 

Data assimilation is a downscaling process for estimating the state of a system by 
synchronizing information from collected coarse mesh measured data and prediction 
derived from numerical forecast model. The classical method of continuous data 
assimilation, see, e.g., [11], is to insert observational measurements directly into 
a model as the latter is being integrated in time. One way to exploit this is to 
insert Fourier low mode observables from a time series into the equation for the 
evolution of the high modes. After a relatively short time t = i, the solution to 
the equation for the high modes is close to the high modes of the exact reference 
solution associated with the observables. At that point the low modes and high 
modes can be combined to form a complete good approximation of the state of the 
system at time t = t, which can then be used as an initial condition for a high 
resolution simulation. This was the approach taken for the 2D Navier-Stokes in 
[6, 7, 21, 22, 26, 29, 30]. Except of the work in [6] (for the 3DVAR Gaussian filter), 
and [5] (which is using the determining parameters nudging approach of this paper 
for data assimilation), the previously mentioned theoretical work assumed that the 
observational measurements are error free. Notably, the authors of [21] present an 
algorithm for data assimilation that uses discrete in space and time measurements. 

In [4] , a new approach was introduced based on an idea from control theory [3] . 
The new approach was motivated by the fact that, unlike the case of finite Fourier 
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modes, measured data represents the values of solution on a discrete set of nodal 
points are not possible to insert directly into the various terms of the evolution 
equations involving spatial derivatives of such solution. In the new algorithm, 
rather than inserting the observational measurements directly into the equations, 
the authors in [4] introduced a feedback control term that forces (nudges) the model 
towards the reference solution corresponding to the observations. This is motivated 
by the fact that instabilities in dissipative evolution equations occur at the large 
spatial scales, hence the need to control, stabilize, or nudge these coarse scales in 
downscaling algorithms. 

The downscaling algorithm for the general setting can be formally described 
as follows: suppose that u{t) represents a solution of some dissipative dynamical 
system governed by an evolution equation of the type 



( 1 . 1 ) 


where the initial data u(0) = uq is missing. Let Ihiuit)) represent an interpolant 
operator based on the observational measurements of this system at a coarse spatial 
resolution of size h, for t G [0,T]. The algorithm proposed in [4] is to construct 
a solution v{t) taking into account the observations that satisfies the evolution 
equations 


— = F{v) - fi{Ih{v) - Ih{u)), 
i;(0) = uo. 


(1.2a) 

(1.2b) 


where /x > 0 is a relaxation (nudging) parameter and vg is taken to be an arbitrary 
initial data. Notice that if system (1.2) is globally well-posed and Ihiv) converge to 
Ih{u) in time, then we recover the reference u{x,t) from the approximate solution 
v{x,t). The goal is to find estimates on /x > 0 and h > 0, in terms of physical 
parameters of the evolution equation (1.1), such that the approximate solution v(t) 
is with increasing accuracy to the reference solution u{t) as more continuous data 
in time is supplied. After some large enough time T > 0, the solution v{T) can 
then be used as an initial condition in system (1.1) to make future predictions of 
the reference solution u{t) for t > T, oi one can continue with (1.2) itself, for as 
long more measurements are provided. 

This algorithm was designed to work for general dissipative dynamical systems 
of the form (1.1). Such systems are known to have global, in time, solutions and 
a finite-dimensional global attractor (for this very reason we are still unable to es¬ 
tablish any results concerning this algorithm for the 3D Navier-Stokes equations; 
nevertheless, the algorithm can still be implemented and tested in practice). More¬ 
over, these systems are also known to have finite set of determining parameters. 
A projection (onto say a finite number of low Fourier modes, or other types of 
interpolant projections based on nodal values and volume elements) is said to be 
determining if, whenever the projection of two trajectories of (1.1) on the global 
attractor approach each other, as t ^ oo, the full trajectories approach each other, 
see, for example, [9, 15-19, 23-25] and references therein. The estimates provided 
on ^ > 0 and h > 0 use the estimates for the global existing solution in the global 
attractor of the system, in terms of the physical parameters, such as Reynolds 
number, in the context of the Navier-Stokes equations, for example. 
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In the context of the incompressible 2D NSE, the authors in [4] studied the 
conditions under which the approximate solution u(t), obtained by this algorithm 
of data assimilation, converges to the reference solution u(t) over time. An extension 
of this approach to the case when the observational data contains stochastic noise 
was analyzed in [5]. Under these assumptions, the data assimilation algorithm 
consists of a system of stochastically forced Navier-Stokes equations. The main 
result established in [5] gives resolution conditions which guarantee that the limit 
supremum, as the time tends to infinity, of the expected value of the L^-norm of 
the difference between the approximating solution and the actual reference solution; 
i.e. the error, is bounded by an estimate involving the variance of the noise in the 
measurements and the spatial resolution of the collected data, h. 

In addition to understanding how the noise in the data affects the accuracy of 
the prediction, another major problem in data assimilation is that typically the 
dimension of the observation vector is less than the dimension of the model’s state 
vector. For example, in order to improve hurricane prediction and typhoon fore¬ 
casts it is important to get the water vapor observations into the system model. 
Water vapor is an important factor in the genesis of tropical cyclones [2] , but accu¬ 
rate and timely measurements still remain a challenge with our current observing 
technology (some recent technological advances in the design of atmospheric water 
vapor observing technology like the GPS Radio Occultation are currently underway 
[36]). In the mean time, it is important to analyze the validity and success of a data 
assimilation algorithm when some state variable observations are not available as 
an input on the numerical forecast model. 

Our main idea in this work stems from the work in [13] which proposes a con¬ 
tinuous data assimilation scheme to the two-dimensional incompressible Benard 
convection problem. The 2D Benard convection problem is given by 


— — vAu + (u ■ V)u + Vp = 002, 

do 

— — kAO -I- (m • V)0 — M • 02 = 0, 
V ■ M = 0, 


(1.3a) 

(1.3b) 

(1.3c) 


u{0,x,y) = uo{x,y), 0(0, x, y) = 0o(x, y), (1.3d) 


with appropriate boundary conditions. Here u > 0 is the fluid viscosity, k > 0 
is the diffusion coefficient, 02 = (0,1) is the second standard basis vector in 
The unknowns are the fluid velocity u{t,x,y) = {ui{t,x,y),U 2 {t,x,y)), the fluid 
pressure p{t,x,y), and the scalar function 0(t,x,?/), which may be the fluctuation 
of the density of the fluid or the temperature of the fluid. The authors in [13] 
proposed an algorithm for the construction of U{t) and ? 7 (t), that approximates the 
velocity u and the temperature fluctuations 0, respectively, from the observational 
measurements Ih{u{t)) of the two components of the velocity field (but without 
needing the measurements Ih{d{t)) for the temperature fluctuations), for t S [0,T] 
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whose evolution is given by 
dU 

— - uAU + {U . V)t/ + VP = 7762 - fiihiU) - Ih{u)), 
ot 

dn 

— - kAtj + {U ■ V)77 - [/ • 62 = 0, 
ot 

V • C/ = 0, 

U{0,x,y) = Uo{x,y), r]{0,x,y) = T]o{x,y), 


(1.4a) 

(1.4b) 

(1.4c) 

(1.4d) 


with the corresponding appropriate boundary conditions; where /i is a positive re¬ 
laxation (nudging) parameter, which relaxes the coarse spatial scales of U toward 
the observed data, P is the approximate pressure, and Uq, rjo are taken to be arbi¬ 
trary. Notice that this algorithm is different from the general algorithm presented 
in [4]. The algorithm in (1.4) construct the approximate solutions for the velocity u 
and temperature fluctuations 6 using only the observational data, Ih{u), of the ve¬ 
locity and without any measurements for the temperature (or density) fluctuations. 
Ideally, one would like to design an algorithm based on temperature measurements 
only, but an algorithm is still out of reach for (1.3). 

In this work, we introduce an abridged dynamic continuous data assimilation for 
the 2D NSE inspired by the recent algorithms introduced in [4, 13]. We establish 
convergence results for the improved algorithm where the observational data needed 
to be measured and inserted into the model equation is reduced or subsampled. Here 
we handle the idealized case where the measured observational data is assumed to 
be free of noise and that the model parameters are exact. To be more precise, to 
review some relevant literature and to set some notation, we start by recalling the 
2D NSE which can be written as 

- vAui + uidxUi + U2dyUi + dx_ 

— 1'Au 2 + UidxU2 -f U2dyU2 + dy 


/l, 

(1.5a) 

/2, 

(1.5b) 

0, 

(1.5c) 

'A{x,y), 

(1.5d) 


where {ui{t,x,y),U 2 {t,x^y)) is the velocity of the fluid at time t and position 
(x, y) € fl, V > 0 represents the kinematic viscosity, pit, x, y) is the pressure and 
{fi{x,y), f 2 (x,y)) is the body force applied to the fluid, which we assume to be 
time independent. We will consider system (1.5) in a physical domain fl, with 
either no-slip boundary conditions or periodic boundary conditions. In the case 
of no-slip Dirichlet boundary conditions we take it = 0 on dfl. The domain D is 
an open, bounded and connected domain of with boundary. In the case of 
periodic boundary conditions we require it, p and / to be L-periodic, in both x 
and y directions, with zero spatial averages over the fundamental periodic domain 

n = [o,t] 2 . 

Inspired by [13], our proposed downscaling algorithm for the construction of ap¬ 
proximate solution, U(t,x,y) from the observational coarse measurements of the 
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one component of the velocity, e.g., the second component, Ih{u 2 {t)), for the refer¬ 
ence solution u(t,x,y), for t € [0,T] is given by 
f)U 

- 12 AU 1 + Uid^Ui + U 2 dyUi + d^P = fi, (1.6a) 

311 

- vAU2 + Uid^U2 + U2dyU2 + dyP = f2- m(4(C/2) - Ih{u2)), ( 1 . 6 b) 

d,Ui+dyU 2 = 0, (1.6c) 

Ui{0,x,y) = U^{x,y), U2i0,x,y) = U^{x,y). (1.6d) 

Here, /r is a positive nudging parameter, which relaxes (nudges) the coarse spatial 
scales of U 2 toward those of the observed data Ih{u 2 ), P is the approximate pressure. 
A choice for and C/^ i® arbitrary, and can be simply taken to be t/° = 0 and 

1 / 2 = 0 . If we knew and U 2 , then we could take C/° = and t/® = ^2 

the solution (Ui,U 2 ) will be identically (ui,U 2 ), by the uniqueness of solutions 
of system (1.6). But, and u® not available, which is the main reason for 
introducing data assimilation algorithms. We note that this algorithm requires 
observational measurements of only one component of the velocity vector field, 
horizontal Ih{ui{t)) or vertical Ih{u 2 {t)). Here, the observational measurements 
Ih(u 2 {t)) were chosen as an example. 

We will consider interpolant observables given by linear interpolant operators 
Ih '■ H^{n) —> L^(r 2 ), that approximate identity and satisfy the approximation 
property 

y - Ih{/>)\\L^(^n) < 7o/i||¥’llffi(n) > (1-7) 

for some positive constant 70 and for every ip in the Sobolev space We 

also consider a second type of interpolant observables given by linear interpolant 
operators =>■ L^(n), that satisfy the approximation property 

y - 4 (<p)llL"(n) < 71^ ll<Pllffi(n) + 72^^ ll¥^lliL 2 (n), ( 1 . 8 ) 

for some positive constants 71 , 72 and for every (p in the Sobolev space H^(fl). One 
example of an interpolant observable that satisfies (1.7) is the orthogonal projection 
onto the low Fourier modes with wave numbers k such that |A:| < 1/h. A more 
physical example are the volume elements that were studied in [24]. An example 
of an interpolant observable that satisfies ( 1 . 8 ) is given by the measurements at 
a discrete set of nodal points in H (see Appendix A in [4]). We will call the 
interpolants that satisfy (1.7) and (1.8) of type I and type H, respectively. 

We provide explicit estimates on the spatial resolution h of the observational 
measurements and the relaxation (nudging) parameter in terms of physical pa¬ 
rameters, that are needed in order for the proposed downscaling algorithm to re¬ 
cover the reference resolution. While the typical scenario in data assimilation is to 
choose p depending on /i, in our convergence analysis we choose our parameters y 
and h to depend on physical parameters. More explicitly, we choose y to depend on 
the bounds of the solution on the global attractor of the system and then choose h 
to depend on y and the physical parameters. The philosophy here is that in order 
to prove the convergence theorems, we need to have a complete resolution of the 
flow, so h has to depend on the physical parameters a.k.a the Grashof (Reynolds) 
number. Numerical simulations in [20] (see also [21]) have shown that, in the ab¬ 
sence of measurements errors, the continuous data assimilation algorithm ( 1 . 2 ) for 
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the 2D Navier-Stokes equations performs much better than analytical estimates in 
[4] would suggest. This was also noted in a different context in [29] and [30]. It is 
likely that the data assimilation algorithm studied in this paper will also perform 
much better than suggested by the analytical results, i.e. under more relaxed con¬ 
ditions than those assumed in the rigorous estimates. This is a subject of future 
work. 

We can extend the corresponding convergence analysis for the 2D Benard equa¬ 
tion, where the approximate solutions constructed using observations in only one 
component of the two-dimensional velocity field and without any measurements on 
the temperature, converge in time to the reference solution of the 2D Benard sys¬ 
tem. This will be a progression of a recent result in [13] where convergence results 
were established, given that observations are known at discrete points on all of 
the components of the velocity field and without any measurements of the tem¬ 
perature. The proposed data assimilation algorithm can also be applied to several 
three-dimensional subgrid scale turbulence models. In [1], it was shown that ap¬ 
proximate solutions constructed using observations on all three components of the 
unfiltered velocity field converge in time to the reference solution of the 3D NS-a 
model. Morever, in [28] a similar data assimilation algorithm was introduced for 
the 3D Brinkman-Forchheimer-Darcy model for flow in porous media. We give a 
progression to this scheme and propose that one can show a sharpened results that 
the approximate solutions constructed using observations in only any two compo¬ 
nents and without any measurements on the third component of the velocity field 
converge in time to the reference solution for this model. The analysis for both 
applications are work in progress. Notably, a similar algorithm for stochastically 
noisy data is at hand combining ideas from the present work and [5]. 

In the section 2, we lay out the functional setting commonly used in the math¬ 
ematical study of the Navier-Stokes equations. The main results are in sections 3 
and 4. We find estimates on the adequate resolution in the observational data h 
and the relaxation parameter /i for observational measurements that satisfy (1.7) 
and (1.8), separately, in the cases of Dirichlet and periodic boundary conditions. 
We prove the well-posedness of system (1.6) as well as the convergence (at an ex¬ 
ponential rate) of the approximate solution v{t) of (1.6) to the reference solution 
u{t) to the 2D Navier-Stokes equations (1.5). 

2. Preliminaries 

For the sake of completeness, this section presents some preliminary material 
and notation commonly used in the mathematical study of hydrodynamics models, 
in particular in the study of the Navier-Stokes equations (NSE) and the Euler 
equations. For more detailed discussion on these topics, we refer the reader to, e.g., 
[10], [31], [33] and [34]. 

In the two-dimensional case with no-slip Dirichlet boundary conditions, let D be 
an open, bounded and connected domain with boundary. We define V to be the 
set of divergence free and compactly supported C°° vector fields from 11 C t 
In the case of periodic boundary conditions, let ft = [0,L]^ for some fixed L > 0, 
we define V to be the set of all L-periodic trigonometric polynomials from to 

that are divergence free and have zero averages. We denote by LP{Q), 1E®’^(0), 
and = VF®’^(0) to be the usual Sobolev spaces in two-dimensions. We will 
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denote by H and V the closure of V in and respectively. We also 

denote by ffo(fi) the set of functions with zero traces at the boundary dil. 

In the periodic case we denote by the subspace of H^{H) of functions that 

are periodic with zero average. 

We define the inner products on L^{H) and H by 

( m , w) = / UiWi dxdy, 

and the associated norm ||u IIl^(O) = Notice that V, and 

are Hilbert spaces with the inner product 



with the associated = ((m, Note that ((•,•)) defines 

a norm due to the Poincare inequality (2.1), below. 

Remark 2.1. We will use these notations indiscriminately for both scalars and 
vectors, which should not be a source of confusion. 

We denote the dual of H by F . Define the Leray projector Pa as the orthogonal 
projection from L^(f2) onto H, and define the Stokes operator A : V —> V is given 

by 

Au = —Pa-Au, 

with domain T>(A) = V D The linear operator A is self-adjoint and pos¬ 

itive definite with compact inverse A~^ : H ^ H. Thus, there exists a com¬ 
plete orthonormal set of eigenfunctions Wi in H such that Awi = XiWi where 
0 < Ai < A 2 < ... < Ai < Ai+i < ... for i G N. 

We also recall the Poincare inequalities: 


(1) For all G y: 





(2.1a) 

(2) for all ip G 2A(a1): 





(2.1b) 


where Ai is the smallest eigenvalue of the operator A in two-dimensions, subject to 
the relevant boundary conditions. 

Let y be a Banach space. We denote by LP{[0,T]-,Y) the space of (Bochner) 
measurable functions t >->• ic(<), where w{t) G F, for a.e. t G [0,T], such that the 
integral \\w(t)\\Y dt is finite. 

Hereafter, c, Cl and ct will denote universal dimensionless positive constants. 
Our estimates for the nonlinear terms will involve the Ladyzhenskaya’s inequality 
in two-dimensions for an integrable function G F: 

( 2 - 2 ) 

Remark 2.2. We note that the Poincare inequality (2.1a) and the Ladyzenskaya 
inequality (2.2) hold for Lp G Hq( 0) or G iL^(D). Similarly the Poincare in¬ 
equality (2.1b) holds for ip G iL^(D) n iLo(^) or G iL^(D) n iL^(D). We will use 
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these versions of the Poincare inequality and the Ladyzenskaya inequality for the 
components of the velocity vector field. 


Also, we will use the following logarithmic estimates for the nonlinear term in 
two-dimensions. These estimates are the analogue of the logarithmic estimates 
proved in [35], for the advection nonlinear term of the Navier-Stokes equations; 
and they can be proved following the same steps of the proof in [35] . 

(1) For every u,v,w G (or H^{n)), with w ^ 0, we have 


u diV w dxdy 


< 


CT ||Vm|]^2(j.j) ]]Vn||^2(j.j) IkllL^(n) fl + log (^rjT, 


1/2 




(2.3a) 


L-‘(n), 


(2) For every u G ^^o(f2) and v,w € H'^(Vl) n i?o(f2) (or u G H^{n) and 
v,w € nH^{n)), with n 7^ 0, we have 


udiV djjW dxdy 


< 


CT IIVn||i2(s.2) IIVn||i2(f2) ||An;||^2(j.;) ("l -b log f 


1/2 


. (2.3b) 


(3) For every v G i^o(^2) and u,w G fl Ho(^) (or u G and 

v,w G ndl^p)), with u 7^ 0, we have 


u diV djjW dxdy 


< 


ct ||Vm||^2^q^ ll^^llL2(n) ll^^llL^(n) 


-b log 


/ ||Au||^2^j.2^ \\ 

\Ai^ II^^IIl2(o) / / 


1/2 


(2.3c) 


where di is interchangeable with dx or dy, and djj is interchangeable with dxx 
or dyy. We note that the logarithmic estimate (2.3c) can follow by an argument 
using the Brezis-Gallouet logarithmic inequality [8]. In [35], it was proven using a 
different approach. On the other hand, the logarithmic estimate (2.3b) does not 
follow as a consequence of the Brezis-Gallouet inequality (see [35] for the proof). 

We recall that in two dimensions and in the case of periodic boundary conditions 
the nonlinearity also satisfies 


((u-V)u,Au) = 0, (2.4a) 

for each u G T>(A) and consequently 

((u ■ V)w, Aw) + ((w ■ V)u, Aw) + ((w ■ V)w, Au) = 0, (2.4b) 


for each u and w G 'D{A). 

Furthermore, inequality (1.7) implies that 

\\w - 4(w)|]i2(j.j) < clh^ ||Vw||^2(fj) , 


(2.5) 
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for every w G V, i/g(r2), or where co = 70 , and respectively, (1.8) implies 

that 


Ik - 4(w)llk(n) < ||Vw||^2(j. 2) + 7Co^'‘ I|Aw||^2 


(n) 


( 2 . 6 ) 


for every w S V^A), i?o(n) D or H^{D,) n iJ^(n), for some cq > 0 that 

depends only on 70 , 71 and 72 . We note that in the case of periodic boundary 
conditions, we demand that the spatial average of Ih(w) to be zero, for every w in 
the relevant domain of Ih (c.f. [4]). This is to guarantee that the spatial average 
of the solution U of (1.6) is preserved, and hence can be chosen to be zero. 

We will use the following elementary inequality proved in [4]. 

Lemma 2.3. Let (j){r) = r — 7(1 + log(r)), where 7 > 0. Then 

min{^(r) : t > 1 } > — 7 log( 7 ). 

We will also apply the following uniform Gronwall’s inequality proved in [24]. 

Lemma 2.4. [Uniform Gronwall’s inequalityj Let t > 0 be arbitrary but fixed. 
Suppose that Y(t) is an absolutely continuous function which is locally integrable 
and that it satisfies the following: 


^+mY<o, 


( 0 ,oo), 


and 


rt+T 


lim inf 

t—¥00 


P{s)ds > 7 , 


ft + T 


lim sup / 
t —>00 Jt 


for some 7 > 0, where fd = max{—/3,0}. Then, Y(t) 
as t —>■ 00 . 


/3 (s) ds < 00 , (2.7) 

0 at an exponential rate, 


Here we denote by G the Grashof number in two-dimensions 

[fTxf II/IIl"(0 ) • 


G = 


( 2 . 8 ) 


We recall that the 2D NSE (1.5) are well-posed and posses a finite-dimensional 
global attractor when / is time-independent, see e.g., [10], [31], [32]. Next, we give 
bounds on solutions u of (1.5) that we will use later in our analysis. These bounds 
are proved in the references listed above. The estimate (2.10b) is proved in [12] 
and estimate (2.9d) will be proved in Appendix A. 

Proposition 2.5. Let r > 0 6 e arbitrary, and let G be the Grashof number given 
in (2.8). Suppose that u is a solution of (1.5) subject to no-slip Drichlet boundary 
conditions, then there exists a time to > 0 such that for all t > to we have 

ft+T 


\Ht)\\]^ < 2kG^ ^ 11 Vm(s)|]^2(j. 2) ^ 2(1 -f tuAi)z/G2 

l|Vu(t)|]^2(j.2) < ckAlGk‘^^ 

J |lAit(s)11^2(fj) ds < (ce*^ -I-tuAi)uAiG^, 


(2.9a) 

(2.9b) 

(2.9c) 


\\^u{t)\\l2^^^ < di^^XlG^ (1 + (1 + Gk^")(l + + Gk*^")) , (2.9d) 


and 
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for some positive non-dimensional constant c. In the case of periodic boundary 
conditions we have 

^ ||Au(s)||^ 2 (j^) c^s<2(l + T^ 2 Al)l/AlG^ (2.10a) 

and 

ll^wWllL2(n) ^ ci/2 A?(1 + G)4. (2.10b) 

Remark 2.6. In this work, we will assume that the reference solution of the Navier- 
Stokes equations, that we are trying to approximate, has evolved enough in time 
to satisfy the estimates provided in the above proposition. That is, we will assume 
that the solution satisfies these estimates at t = 0. For that reason, we will take, 
without loss of generality, to = 0. 


3. Convergence analysis with observable data of type I 

Following similar techniques introduced for the two-dimensional Navier-Stokes 
equations (see, e.g., [10, 31, 34]), we can prove the global well-posedness of system 
(1.6) as stated below, when the observable data satisfy (1.7). For more details, see 

[4]. 

Theorem 3.1. [Well-posedness of solutions] Suppose Ih satisfy (1.7) with p > 0 
and h > 0 are chosen such that pc^h'^ < v , where cg is the constant in (1.7). If the 
initial data Uq £ V, then the continuous data assimilation system (1.6), subject to 
Dirichlet or periodic boundary conditions, possess a unique global strong solution 
U[t,x,y) = (Ui(t,x,y),U 2 {t,x,y)) that satisfies 

U €Ci[0,T]-V)nL^{[0,T]-V(A)), and ^ e L^i[0,T]; H). 

Moreover, the solution U{t,x,y) depends continuously on the initial data Uq. 

We will now state and prove that under certain conditions on p and h, the 
solution {Ui,U 2 ) of the data assimilation system ( 1 . 6 ) converges to the solution 
{ui,U 2 ) of the two-dimensional Navier-Stokes equations (1.5), subject to periodic 
or Dirichlet boundary conditions, respectively, as t —>• oo, when the observable 
operators satisfy (1.7). 

Theorem 3.2. Suppose that Ih satisfy the approximation property (1.7) and u(t, x, y) 
= {ui(t,x,y),U 2 {t,x,y)) is a strong solution in the global attractor of (1.5) sub¬ 
ject to Dirichlet boundary conditions. Let U{t,x,y) = {Ui{t,x,y),U 2 (t,x,y)) be 
a strong solution of (1.6), subject to Dirichlet boundary conditions. If p > 0 is 
chosen large enough such that 

p > 2cvXi{l + log(G) -h G^)G^ 

and h > 0 is chosen small enough such that pcgh^ < v, then \\u{t) — t 7 (t)||^ 2 ^Q^ 
at an exponential rate, as t ^ oo. 


(3.1) 
^ 0 
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Proof. Define u = u — U and p = p — P. Then ui and U 2 satisfy the equations 
du 

— -uAui + UidxUi + U 2 dyUi + uidxUi + U 2 dyUi + dxP = 0, (3.2a) 

C/ L 

du 

— -uAu2 + UidxU2 + U2dyU2 + UidxU2 + U2dyU2 + dyP = -plh{u2), (3.2b) 

dxUi + dyU 2 = 0. (3.2c) 

Since mi, U 2 ^ ^ and ^ are bounded in L^([0, T]; iJ), we can take the L^(ri) 
inner product of (3.2a) and (3.2b) with ui and U 2 , respectively. We obtain, using 
the divergence free condition (1.6c), integration by parts and using the relevant 
boundary conditions, that 

ll^illz,"(n) +^"l|Vui||i 2 (j.j) < \Jia\ + \Jib\ - {dxP,Ul), 

< \ J2a\ + \J2b\ - {dyP,U2) - p{Ih{u2),U2), 

where 


Jla ■= {uidxUi.Ui), Jib := {u 2 dyUi,Ui), 

J2a ■= {uidxU2,U2), J26 := (U2l9yU2, U2) • 

To estimate the nonlinear terms we proceed as follows: using integration by parts 
twice, we have 

Jla = {dxUi, (ui)^) = -2 [uiudxUi) 

= 2 {uiUi,dyU2) 

= -2 {uidyUi,U2) - 2 {uidyUi,U2) 

= : —2(Jiai) — 2( Jia2), 


where we used above the divergence free condition (3.2c) and the relevant boundary 
conditions. By the logarithmic estimate (2.3a) and Young’s inequality, we have 

\Jlal\ '■= I {uidyUl,U2) I 

" " iivu 2|L2(^) 


< CT ||Vui||i2(f2) l|VMi||i2(j.;) \\u2\\L^(^n) 1 + log ^ 


yr\\u2h^ 


Lqn), 


U 2 , C 2 K . , f Il^^2|li2(j2) 

- ^ + - ll^«illL’^(n) I 1 + log ' 


^ 


(3.3) 


and 

|<^la 2 | := I {uidyUi,U 2 ) I 




||Vu 


2 ry 


L^{n) 


L-‘{Q.) 


A^"\\n2\ 


L^(n). 


II'“2 |Il 2 (j.j) . 

(3.4) 


1 + log 
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Thus, 


\Jla\ < ll^“lllL^(n) + ll^''^lllL^(n) I 1 + log I 1/2 


IIVU2II 


L^(n) 


K \\U2\ 


L^(n), 


(3.5) 


A similar argument as in (3.3) and (3.4) yields 
\Jlb\ ■= \{U2dyUl,Ul)\ 

V 2 , C 2 A , , i 11 ^^ 211 ^ 2 ( 5 . 2 ) 

< ^ I|Vmi|Il 2 (s 2 ) + - ||Vui||i^ 2 (s 2) 1 + log ' 


32 

and 

\J2a\ '■= \{uidxU2,U2)\ 


- ^ + I Il^«2|li2(5.2) ( 1 + log 


A^^\\n2\ 


L-'(n). 


IIVU2II 


L=(n) 


AA\\n2\\ 


L^(n), 


Il^2|li2(j.2) . 

(3.6) 


ll'“2|li2(Q( . 

(3.7) 


By Holder’s inequality, Ladyzhenskaya’s inequality (2.2), we have 
1^261 = \{u2dyU2,U2)\ < 11^211^4(51) ||9yU211^2(51) 

< Cl ||u2||L2(n) I|Vc(2||l 2(51) ll^yC(2||i^2(5i) 

“ ^ ll^'“2llL=(n) + - 11 '“2 11^2(51) ||u2|li2 


(n) 


(3.8) 


Thanks to the assumption fj,CQh‘^ < v and Young’s inequality, 
-fj.{lhiu2), U2) = -^i{Ih{u2) - U2, U2) - M l|w 2 |lL 2 (n) 

< mII4(w2) -M 2|lL=m) l|w2|lL2(n) -Ml|u2|li2 


^(n) 


< /iCoh ||M2|lL2(n) l|VM2|lL2(n) - M I|u 2 |Il 2 


^(n) 


|2 _ /^ 


< ||V{t2||ji255ii ^ l|w2|lL2(n) 


< 2 Il^“2|li2(5i) - I 117121112(51) . 


Also, note that 


(3.9) 

(3.10) 


{ dxP , Ml) + { dyP , U2 ) = 0, 

thanks to the divergence free condition (3.2c), integration by parts and using the 
relevant boundary conditions. It follows from the estimates (3.5)-(3.10) that 


111^11^(0 ) + ^ I|Vm||12(5i) < 


^ IIVm||12(51) 1^1 + log )) “ ^) II“2 |I^(q) , 


iiL4(n), 
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or 




4 

c ,,2 


- 1^1 +log y— 


2 

L=(n), 


M l|U2||i2(j^) 


(3.11) 


Next, we will use a similar argument as in [4] and [14] to prove that ||w(t)U2(o) 0; 

exponentially, as t —>■ oo. 

Using the Poincare inequality (2.1), we may rewrite (3.11) as 


“ llM|lL2(n) + 


vX 


dt 


1 ll~l|2 , 


^Ai I|V'U2 |Il2(s.j) 

4 Ai II“2 |Ii,2(j.j) _ 


^ l|VM||i 2 (j.;) M+log 


^2|lL^(n) 
2 


211^2(0) 


Ai 


2 

L^(n), 


-/i ||u 2 |li 2 


L^{Q) 


or 


vXi 


vXi 


dt 


4fi\ 


|M|^^2,n^ + — l|M|lL2(n) + — (^'(>(»'(l)) + ^ j 11^2 |li2(s.j) < 0, 
where we denoted by 

^(r(t)) := r{t) - 7 (t)(l + log(r(t)), 

l|V'“2|li2(j.2) c 2 

"(0 := ' „2 -. 7(U := yru • 


•^1 II«2|Il2(5.2) 

Now, Lemma 2.3 implies that 


v'^Xi 


(3.12) 

(3.13) 


(j){r{t)) > -7log(7) = l|Vw||i 2 (s 2 ) log IIVw||i 2 (fj)^ . 

Setting 

Pit) l|VM||i 2 (fj) log l|Vu||^ 2 (s.j)^ , 

we have 

^ I|ll|li2(n) + (l|llillL=(n) + II'“2 |Il2(q)^ + Pit) 1111211^2(5.2) ^ 0. (3-14) 

We may conclude that 

^ l|w|li=(n) + min{^. Pit)} ||w||i2(s.2) < 0. 

Taking r = (j^Ai)”^ in Proposition 2.5, using (2.9a) and (2.9c), and due to Remark 
2.6, we conclude that 

^ ^ l|Vn(s)||^ 2 ( 52 ) log l|VM(s)||i 2 ( 5 ^(^ ds 

< -(1+log(G) TCJ'^) [ ||Vu(s)||i 2 (s. 2 ) ds 

^ Jt 

< c(l + log(G) + G^)G^ 


(3.15) 
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for all t > 0. Therefore, the assumption (3.1) implies that 

Pt + T 


lim inf 

t—>oo 


/3(s) ds > 




2z/Ai 


> 0 , 


and 

3/1 

limsup / f3{s)ds < -—— < oo. 

t-J-oo Jt 2v\i 

Define P{t) := min{^^, /3(t)}, then P{t) satisfies the condition (2.7). 
uniform Gronwall’s lemma (2.4), we obtain 

at an exponential rate as t —^ oo. 


(3.16) 


(3.17) 
By the 


□ 


Remark 3.3. Theorem (3.2) holds in the case of periodic boundary conditions. 
Using (2.10a) in Proposition 2.5, the estimate (4.31) can be improved to 

^■s^c(l + log(G))G^ (3.18) 

thus the lower bound on /i (3.1) can be improved to 

/I > 2ciiAi(l + log(G))G^. (3.19) 


4. Convergence analysis with observable data of type 11 

Next, we will prove that under certain conditions on /i, the approximate solution 
(Gi, U 2 ) of the data assimilation system (1.6) exists globally in time and converges 
to the solution (111,112) of the 2D Navier-Stokes equations (1.5), subject to periodic 
or Dirichlet boundary conditions, respectively, as t —>■ 00, when the observable 
operators satisfy (1.8). 

4.1. Dirichlet Boundary Conditions. The existence and uniqueness of strong 
solutions of the data assimilation system (1.6) with observables that satisfy (1.8), 
subject to Dirichlet boundary conditions, does not follow immediately as in the 
case of the 2D Navier-Stokes equations. Extra conditions on the nudging constant 
/I > 0 are required. Next, we will prove the global existence of the strong solution 
U{t) of (1.6) and simultaneously show that it converges, in time, to the reference 
solution u{t) of the 2D Navier-Stokes equations (1.5). 

Theorem 4.1. Let Ih satisfy the approximation property (1.8) and u{t,x,y) = 
(ui(t, cc, y), it2(t, a:, t/)) be a strong solution in the global attractor of (1.5) subject to 
Dirichlet boundary conditions. If £ V with 

||VC/°||"2(f^) <5i^"AiG2e^", (4.1) 

where c is the constant in (2.9b), and p > 0 is large enough such that 

p > 2cvXiK\og{K), 


(4.2) 
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where K is defined in (4.18), and h > 0 is chosen small enough such that gic'^h'^ < v. 
Then, there exists a unique strong solution {Ui{t,x,y),U 2 {t,x,y)) of (1.6) that 
satisfies 

U €Ci[0,T];V)nLm0,T]-ViA)), and ^ e L\[0,T]; H), 

such that 

||VC/(t)||i.(^) < Tcu^MG^e^\ 

for all t > 0. Moreover, ||Vu(t) — VC/(i)||^ 2 (-Q^ 0, at an exponential rate, as 

t —>• oo. 


Proof. We will prove some formal apriori estimates that are essential in proving 
the global existence of solutions of system (1.6). These estimates can be justified 
rigorously by using the Galerkin method and the Aubin compactness theorem (see 
e.g. [10]). 

Define u = u—U andp = p—P. Then ui and U 2 satisfy the equations (3.2). Since 
by assumption, ( 111 , 1 / 2 ) is a solution which is contained in the global attractor of 
(1.5), in particular, it satisfies the global estimates in Proposition 2.5, then showing 
the global existence, in time, of the solution {ui,U 2 ) is equivalent to showing the 
global existence, in time, of the solution {Ui, U 2 ) of system (1.6). To be concise here, 
we will show the global existence of the solution u{f) and show that ||^^(t)||i 2 ^Q^ 
decays exponentially, in time, which will prove the convergence of the approximate 
solution Uif) to the exact solution u{t), exponentially in time. 

Since || < ciy^XiG^e^‘', then by the continuity of ||V17(t)||^2(Q^, there 

exists a short time interval [0,T) such that 

mmlfin) < (4.3) 

for all t € [0,T). Assume [0,T) is the maximal finite time interval such that (4.3) 
holds. We will show, by contradiction, that T = 00 . Assume that T < 00 , then it 
is clear that 

limsup ||V17(t)||^2^Q^ = Icu^XiG^e^ , 
t^f- 

otherwise, (4.3) holds beyond T. Taking the inner product of (3.2a) and 

(3.2b) with —Ami and —Am 2 , respectively, we obtain, on the time interval [0,T), 
that 

2 dt ^ ^ \ J^a\ + l^lhl + \ Jlc\ + l^lhl + [dxP, Ami), 

2 fit + V ||Am2||^2(q) < I J2a| + IJ 26 I + I J2c| + \ J2d\ + {dyp, AM 2 ) 

- pL{Ih{u2), -Am2 ), 

where 


Jla '■= [UidxUi, 
Jic := (uidxUi, 
J2a '■= {UldxU2, 
J2c ■= iuidxU2, 


Aul), 

Jib 

Ami), 

Jld 

Am2), 

J2b 

Am2), 

J2d 


{U2dyUi, -Aui), 
(M 29 yMi, -Aui), 
{U2dyU2, -AM2), 
{U2dyU2, -AM2). 
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Thanks to the assumption < v, Young’s inequality and inequality (1-8), 

we have 


-n{lh{u2 ),-Au2) = -/i(//i(u2) - U2,-Au2) - /i ||Vm2||^2 


(O) 


< M||4(m 2) ||AM2|li2(j^) -/l||u2||i2(j^) 


ly 


< — \\Ihiu 2 ) - U2\\l2,^) + 7 \\AU2\\12,^. - ^ ||Vu 2 ||i 2 




,,2 4l4 

^ ||A« 2 |li 2 (o) + 4 


21 / 4 ^ 

-Ai||Vi( 2 |li 2 (f 2 ) 

<^||A«2|li2(^)-|||Vu2|li2(^). (4.4) 

To estimate the nonlinear terms we proceed as follows: using the divergence free 
condition (3.2c) and the logarithmic inequality (2.3b): 

Jia = {Uid^ui,-Aui) = {UidyU 2 ,Aui) 

/ ^ iia« 2 |L 2 (^) 




L^{n) 

2 

L^{Q) 


< CT ||VC/i||^ 2 ( 5 .j) l|Vu 2 |li 2 ( 5 .j) ||Aiii||^ 2 (j.j) I 1 + log I 


Ai^^||Vzi 2 lL 2 (^), 

I|a«2|| 


- Too + Z ll^^illi’^(n) W^Ml^q) 1 + log ^ 


L^(n) 


yr\\vu2h- 


Lqo), 


(4.5) 


By a similar argument, we can show 
J2a = (UldxU2, —AU 2 ) 

<^||Afi 2 |li 2 (^) + ^||VC/i||^ 2 (^)||Vli 2 |li 2 (^) fl + log^ I|A^^ 2 |L 2 (^) 


100 

and 

J2b = {U2dyU2,—Au2) 


\\'^\\Wu2\\l^ 


L^(n). 


(4.6) 


- ^ ||VC/ 2 ||i 2 (^) l|Vii 2 |li 2 (^) (1 + log ( 




(4.7) 


By the logarithmic inequality (2.3c), we have 

Jid = {u2dyUi,-Aui) 

< CT ||Vii 2 ||i 2 (j. 2 ) IIVui 11^2(5.2) l|Aili|li 2 (j.j) 1 + log 


I|a« 2 || 


1/2 


L=(n) 


A}/^llvu 2 iL 2 (^)yy 

l|AM2|li2^Qj 


“ 100 jy II^'“ 2 |Il 2 (q) I 1 + log 


A}/"||Vti 2 lL 2 (^), 


(4.8) 
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Similarly, we have 
J2d = iu2dyU2 ,-Au2) 


< ||A{i2||i2(s^) + - ||VU2|li2(s^) ||VU2 |Il 2 (s^) 1 + log _ 


||Am: 


2 |Il 2 (q) 


Ar I|v« 2|L2(^), 


(4.9) 

Integration by parts and the boundary conditions yield 

Jib — (^U2dyUl^ Aril) = {U2dyUl’) dxx’^l') 4“ (^U2dyUl^ dyyUl') 

— (^dxU2dyUl ^ dxUl'j “t“ (U2dyx‘^l: dx'U'l') (U2^yUl ^ dyyUl'j 

= ■ Jlbl + Jlb2 — Jlb3- 

Using the divergence free condition (3.2c) and the logarithmic inequality (2.3a), we 
get 

Jlbl — i^xU20yUl^ dx^l) = {dyUldxU2T dyU2) 

1/2 


< CT 11 Aril 11 ^2 \\VU 2 \\i^ 2 (q) ||VM2||i2fj.j^ I 1 + log 


l|Afi 2 || 


L^{Q) 


L^{n) II'' 


~ 100 ^ ||VM2||i2(f;) f 1 + log 


Xy"\\Vu2h^^n)J J 

^ ||Au2|li2(s^) 


,Ai^ ||Vm 2||^2 ^j. 2) / / 
(4.10) 

Similarly, the divergence free condition (3.2c) and the logarithmic inequality (2.3b) 
imply 

Jlb2 = {U2dyxUl,dxUl) = {U2dyyU2,dyU2) = {U2dyU2,dyyU2) 

( ( ||Au2|li2.f^. 

<ct||V[/2|| ||Vfi 2 || L^(a) l|A'«2||i2|.Qj I 1 + log I ^^2 


“ 100 ||VM2|li2(j.j) f 1 + log 


Ar^iiv«2iL2(^)yy 

l|Au2|li2(s^) 


,Ai^ I|V'«2||^2 ^qj / y 
(4.11) 

By integration by parts and the divergence free condition (1.6c), we also have 

Jibs — {U2dyUlj dyyUl') = ~(U2^y{^y’^l) ) 

= -^idyU2idyUi)^) = ^{dxUiidyUif) 

— ) — (Ul9y^l^yy‘^2) 

— ~{9yUldyUl^ dyU2) ” (y {t 1 7/2 ) =: ^IbSl '^1632- 
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Following similar argument as above, using the logarithmic inequalities (2.3a) and 
(2.3b), we can show that 

<^1631 ~ (^dyUldyU\^ dyU2') — iJ^yU\dyU\ ^ dyti^^ 


< 


100 


II Awi 11^(0) + - l|Vt/i||^..(o) l|Vu2|li2(^) 1 + log 


I|A«2|| 


L=(n) 


(4.12) 


and 

Jlb32 — (UldyyUij dyU2^ — (UldyU2 j ^yy^l^ 
v ^ ..2 C 


< 


100 


l|A«i||^(0) + ^ l|Vt/i||^(o) l|Vu2||^(o) ( 1 + log ( 


(4.13) 

Integration by parts, the divergence free condition (3.2c) and the boundary con¬ 
ditions imply 

Jlc — ^lll) — 4“ ^yy’^l') 

= -{dxUldxUl,dyU2) - {UldxxUl,dyU2) - {UldxUl,dyyUl) 
= '■ Jlcl + Jlc2 + JlC 3 - 
By the logarithmic inequality (2.3a), we have 
Jlcl — {dx^l^X^lT dyU2^ 

/ / \ \ 1/2 

I|Au 2|| \\ 


< c||Aui||^ 2 (j. 2 ) IIVui 11^2(5.2) \\^ML^(n) U + log TT 72 


- iM ^ ll^«ill^(n) ( 1 + log 


V l|Vu2lL=(n)/, 

IIAU2I 


,A}/"||Vu2|L.(n), 


lL"(n) II V «2 |Il^(0) 

\ \^1 l|VM2||i= 

(4.14) 

Applying the logarithmic inequality (2.3b) and the Poincare inequality yield 
Jlc2 = -{uidxxUl,dyU2) = -{uidyU 2 ,dxxUl) 

/ / \ \ 1/2 

IIAU2II \\ 


< c||Aui||^2(s^) l|Vui||i2(s^) l|Vu2||i2(s^) I 1 + log I ^ 


Il^'ii2|li2(n) 


^ - "2 . oA]^ ^ ..2 ||__ ||2 (ill ( IIAM 2 II 


- 100 + ~~;r~ ii^“2||l2(s^) i i + log 


Al/2||Vu2||i2(f,) 

(4.15) ' 
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The term Jic^ can be estimated by a similar argument as above. Also, using 
integration by parts and the divergence free condition (1.5c), we get 

J2c := {uidxU2, -Au 2) = {uidxU2, -dxxU2) + {uidxU2, -dyyU2) 

= idxUldxU2, dxU2) + {uidxxU2,dxU2) 

T {^dyU\dx'^2i ^y'^2) (u\dxy'^2'^ 

= {dxUldxU2,dxU2) + iuidxxU2,dxU2) 

“t“ (^0y'U\dx^2^ ^y'^2^ 1 ^^^ 2 ) 

= : J2CI + J2c2 + J2c3 + J2c4- 

The above four terms can be also estimated using the logarithmic inequality (2.3a) 
as shown previously. 

We note that 

{dxp, Aui ) + {dyp, Au2) = 0, (4.16) 

due to the divergence free condition (3.2c) and integration by parts. It follows from 
the above estimates that on the time interval [0,T), we have 

Using the Poincare inequality, we conclude that 

^ ^ llVwill^(^) + m ||V{t2||^(n) < 0, (4.17) 

with 

/3(0 = ^[r- 7(i)(l + log(T))] + /i, 

where we denoted by 

l|Aw2|li2^Q^ 

Ai||Vu2||i.(n)’ 

and 

lit) = ^ (llVii(t)ll^(^) + ||VC/(t)||i2(^) + Ar' ||A«(t)||)2,(^^) . 

On the time interval [O,!"), we have || VC/||^ 2 ^j. 2 j < Scu^AiG^e'^^. Thus, by (2.9b) 
and (2.9d), we have 

lit) < (scu^AiG^e®* + cu2Ai (l + (l + G^e^") (l + + G^e^"))) 

< cG^ (1 + (1 + (1 + + G^e^")) =: AT, (4.18) 

for all t € [0,T). By Lemma 2.3, we may conclude that 

Pit) > -uAi ( 7 (t) log (7(i)))) + A 

> p — cvXiK \ogiK). (4.19) 


4(n) 






20 


ASEEL FARHAT, EVELYN LUNASIN, AND EDRISS S. TITI 


Therefore, the assumption (3.1) implies that j3{t) > ci'XiK\og{K) > 0 for all 
t € [OjT). Define ^{t) := min{^, /3(t)} > 0, then we can write (4.17) as 

Using Gronwall’s Lemma, we have 

< l|Vfi(0)||i.(o) (4.20) 

for all t G [0,f). Since ||Vfi(0)||^.(^) < 2 ||Vu(0)||^.(^) + 2 ||Vt/(0)||^.(^), then by 
(2.10b) and (4.1), and due to Remark 2.6, we have 

l|Vfi(t)||i2(^) < ||Vu(0)||i.(fj) < (4.21) 

for all t G [0,T). This implies that \\'VU(t)\f^ 2 f^^-^ < Gcv'^XiG'^, for all t G [0,T). 
This in turn will yield a contradiction since 

Icv^XiG^e^ = limsup ||Vt/(t)||^2(j.jj <%cv^XiG^e^ . 
t^f~ 

This proves that T = oo. Thus, the solution u{t) exists globally in time and it 
satisfies 

||Vfi(t)||^(f,) <4^^.2AlG2e«^ (4.22) 

for alH > 0. Following the techniques that were introduced to prove the existence 
and uniqueness of solutions for the Navier-Stokes equations (see for example, [10], 
[32] and [34]), we can show the existence of the solution u{t) of system (3.2), that 
will inherit the estimate (4.21). Moreover, ||fi(t)||£2^Qj decays exponentially in time 
and it inherits the inequality (4.20). The uniqueness and the well-posedness will 
follow by a similar argument as above. Since, by assumption, u{t) is a strong 
solution in the global attractor of the 2D Navier-Stokes equations, then this proves 
the global existence and the uniqueness of the solution U{t) of system (1.6), that 
satisfies \\'X7U{t)\\'^^2f^Q^ < Icu'^XiG'^e^^, for all t > 0. Moreover, we have 

\\Vu{t)-VU{t)\\ G{q.) 0, 

at an exponential rate as t —^ oo. □ 


4.2. Periodic Boundary Conditions. The existence and uniqueness of strong 
solutions of the data assimilation algorithm (1.6) with observables that satisfy (1.8), 
subject to periodic conditions, as stated below, follows by a similar argument as for 
the two-dimensional Navier-Stokes equations. See [4], [10], [31] and [32] for more 
details. 

Theorem 4.2. Suppose Ih satisfy (1.8) and p, > 0 and h > 0 are chosen such 
that < V, where cq is the constant in (1.7). If the initial data Uq G V, then 

the continuous data assimilation system (1.6), subject to periodic boundary con¬ 
ditions, possess a unique global strong solution U(t,x,y) = (Ui(t,x,y),U 2 {t,x,y)) 
that satisfies 

U €Gi[0,T]-V)nL^{[0,T]-ViA)), and ^ € L\[0,T]-, H). 

Moreover, the solution U{t,x,y) depends continuously on the initial data Uq. 
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Theorem 4.3. Suppose that It satisfy the approximation property (1.8) and u(t, x, y) 
= {ui{t,x,y),U 2 (t,x,y)) is a strong solution in the global attractor of (1.5) subject 
to periodic boundary conditions. LetU(t,x,y) = {Ui(t,x,y),U 2 (t,x,y)) be a strong 
solution of (1.6), subject to periodic boundary conditions. If fi > 0 is chosen large 
enough such that 

g,>2cvXi{G^ + G^), (4.23) 

and h > 0 is chosen small enough such that yc^h^ < v, then, ||VM(i) — VC/(1)||^2 ^q^ 

—>■ 0, at an exponential rate, as t ^ oo. 

Proof. The proof follows by similar means as in [4]. Taking the inner product of 
(3.2a) and (3.2b) with —Aui and —Au 2 , respectively, adding the two equations 
and then using the orthogonality property (2.4b), we get 

^ ll^“lli,= (n) <!((«• Am)| - n{Ih{u2), -Am 2), 

where u is the reference solution of the 2D Navier-Stokes equations (1.5). 

Notice that, thanks to the divergence free condition (3.2c) and the boundary 
conditions, we have 

((m • V)m, Am) = [uidxU, Au) + (u 2 dyU, Au) 

= {uidxUi, Aui) + {uidxU 2 , Au 2 ) + {u 2 dyu, Au) 

= — {uidyU 2 , Aui) + {uidxU2, Au 2 ) + {u 2 dyU, Au) 

= ■ Jla + Jib + J2- 

Using the logarithmic inequality (2.3b) and the Poincare inequality (2.1), it 
follows that 


\Jla\ = \{uidyU 2 ,Aui)\ 

< CT IIVmi 11^2(5.2) l|Vw2||i2(s.2) l|Aui||i2(s.2) fl + log ( 

V V^l l|V«2||i2 


1/2 


L^(n), 


< ||Ami||^ 2 (s.j) ||VM 2||i2(s.j) ||AMi||^2(fj) I 1 + log I „ 

V VA/ l|VM2|li2(j.j)^ 


I|Am: 


2llL2(n) 


1/2 


By Young’s inequality we get 


\Jla\ < ^ ||Ami||^ 2 (s. 2) + — ||AMi||^ 2 (fj) l|VM 2 |li 2 (fj) ( 1 + log 


I|Am 2||) 


Ai iivM2ii^(n);; 

(4.24) 

for some positive dimensionless constant c. Similarly, by (2.3b), we can show that 


\Jlb\ < ^ I|Ami||j^2^j.2^ + ||Au2||^2^j.2^ ||VU2||^2^j.2j 


11 Am: 


2llL’^(n) 


, Al ||VM2 ||j;^2(.qj ^ 


(4.25) 


1 + log 
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and by (2.3c), we have 

IJ 2 I < J ||Vu 2 |lL 2 (fj) ^1 + log 


Ai ||Vm: 


for some positive dimensionless constant c. By (4.4), we also have 
-f^{h{u^),-Au-) < ^ ||Ali-||i2(^) - I ||Vfi5||i2(^ 
Therefore, by the estimates (4.24)-(4.27), we obtain 


(4.26) 

(4.27) 




II A l|2 /^1 , 1 f 

-WMlL^n) 1 + log ^ 


:))-) 




,Ar^||Vfi2|L2(^), 

Using the Poincare inequality (2.1), we may rewrite the above inequality as 


_d 

dt 


,2 t/Ai 

iL^in) + — 


IIV7~I|2 , “^^'1 IIV7~I|2 , ^Al II^“2 |Il2(q) ^ 

+ A ll^'“llL2(n) + 4 , ||_- ||2 11^^211^2(0) < 

4 Ai ||VM2|lL=(n) 


^||A«||i 2 (^) (l + log' 


or 


vAi 


12 A 1 


^Ai ||Vfi2|li2(j.2), 




-M l|VU 2 |li 2 (s. 2 ) , 


(4.28) 


dt 


df 2 ) 


17 I|V«IIl 2 ( 0 ) + \\^u\\h(n) + - 7 ^ ( + if ) l|V« 2 ||i 2 (f,) < 0, (4.29) 


j^Ai 


df 2 ) 


where we denoted by 


^(r(t)) := r{t) - 7 (t)(l + log(r(t)). 


II Aw2 11^2^5.2) ^ ^ C 2 

:= ' „2 -> 7(i) := Tm ll^“llL=(n) • 


Ai ||V'U 2 || 2 ^ 2 |.j.j^ 

Now, Lemma 2.3 implies that 


ly^Xi 


Hrit)) > -7log(7) = ll^w|lL"(n) log ll^^l 


2 

L^(n) 


Setting 


Pit) :=Ai-^l|Au||i 2 (j. 2 )log 


^l|A«f 


L^{n) 


we have 


dt 

We may conclude that 


iiv«iii2(^) + ^ \\vu\\i.^^^+m iivu2iii2(^) < 0. 


^ IIVu||i 2 (fj) + min{^, /3{t)} IIVu||^ 2 (j. 2 ) < 0. 


(4.30) 
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Taking t = (uAi) ^ in Proposition 2.5, using (2.10a) and (2.10b), and due to 
Remark 2.6, we conclude that 

||Au(s)||^ 2 (f 2 )log(^^^||Au(s)||i 2 (s^)^ ds 

<^(1 + G)^ ||Au(s)||^ 2 (s^) ds 

< c{l + G)G^ =c{G'^+ G^), (4.31) 

for all t > 0. Therefore, the assumption (4.23) implies that 

liminf / f3{s)ds > ^ > 0, and limsup / /3(s)ds < < oo. (4.32) 

Jt 2 i—)-oo Jt 2 

Define ,3(t) := min{^^, /3(t)}, then /3(t) satisfies the condition (2.7). By the 
uniform Gronwall’s lemma (2.4), we obtain 

l|Vu||^(n) = l|Vu(t)-VC/(t)||i2(^)^0, 

at an exponential rate as t —^ oo. □ 


5. Appendix A 


In this appendix, we will prove estimate (2.9d) following the calculations in [31]. 
The proof is formal. It can be done rigorously using the Galerkin approximation 
and then passing to the limit. 

It follows from the 2D Navier-Stokes equations (1.5) that 

Ou 

< ^ ll^w|lL2(a) + II(m • V)u||i 2 (s^) + . 

L=(n) 


dt 


Using the Ladyzhenskaya inequality (2.2), we can show that if u G D{A) then 


(u ■ V)u||^2,s^^ < Cl ||Vu|| 


2(n) ^ -i ii“iil2(o) 
Using Young’s inequality, we have 


L^(n) 


||Au 


, 1/2 

'G{q,) • 


(5.1) 


du 

dt 


3u 


L^(n) 




The estimates (2.9a), (2.9b) and the definition of the Grashof number (2.8) imply 
that there exists a time tp > 0 such that for t >to 


du 


dt 


L^(Q) 


< cv^ ||Am||^2^j.2^ + cG^ [v^XiG^e^ ^ + cu^A^G^, (5.3) 


for some positive non-dimensional constant c. 

Moreover, integrating inequality (5.2) on the time interval [t+ t] for some t > 0 
and using the definition of the Grashof number (2.8) yield 


rt+T 


< c \ u 


du 

'dt 




ds 


G(q.) 


||Au(s)||i;_2(q) ds-I-- ^ ||Vu||j;^2(q) ds-I-ru'^AiG^^ 


for all t > to- 
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Taking t = (i/Ai) ^ in Proposition 2.5 and using the estimates (2.9a), (2.9b) 
and (2.9c), we get 


t + T 


du 


ds 


L^(n) 


( j-t+T pt + T 

\\Au{s)\\l2^^^ ds + iyG"^ ds + ri/'^A^G 

< c |lA«(s)||i2(f,) ds + l|Vu||i2(n) ds + ti.'^XIG 

< c ^v^AiG^(e‘^ + rvAi) + v^Ai(l + ri/Ai)G®e‘^ + ri^^A^G^^ 

(< 


<cv^AiG^ (l + e*" +GV 


(5.4) 


We now differentiate the equations (1.5) w.r.t. to t and take the inner product 
with ^ to obtain 

l_d 
2 dt 


du 

2 

+ 

du 

2 

< 


du\ 

U, 

dt 

L^{Q) 

dt 

L^{n) 

\\dt J 

’ dt) 


< Cl IIVrt||^2(-Qj 
2 


du 



'm 

L^{Q) 

dt 


V 

< - 
- 2 


^du 


L^(a) 




llL=(n) 


L=(n) 

du 


dt 


L=(n) 


where we used the Ladyzhenskaya inequality (2.2) and Young’s inequality in the 
last two steps. Integrating the inequality on the time interval [s,t + r] and using 
(2.9b) we get 

2 


du . . 

2 

du, , 

»<*+"> 

< 

G{n) 



2 „t+T 

2„G* / 


+ a/AiG^e‘- 


L^{n) 


du 

'dt 


(0 


dL 


L^(n) 


for all to < t < s < t + T. Integrating once again w.r.t. to s on the time interval 
[t, t + t] and using (5.4) yield 


du 

'm 


(t + r) 


< (1 + cvXitG^c^ 


L^(Q) 


i+T 


< ci/^AiG^ (1 + cvXitG^c^ ) (1 + + G^e 


du 

'm 


(s) 


ds 


L’^(n) 


'(4„G 




Since r = (j^Ai) then 


du 

'm 


{t + T) 


2 0“ 


< cv^XiG^ 1 + G^e 




')■ 


L^{n) 

for all t > to- We can redefine to to be large enough such that for all t > to 
2 


du 


dt 


it) 


< ci^^XlG^ ( 1 + G^e^') (l + G^e^') . 

The 2D Navier-Stokes equations (1.5) imply that 

du 


(5.5) 


'IIA^II 




dt 


L^+l) 


+ ||(m • V)u||j^2j-q^ + ||/||^2(f2) ■ 
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This, using the estimate (5.5), the nonlinearity estimate (5.1), the estimates (2.9a) 
and (2.9b) and the definition of the Grashof number (2.8), we can conclude that 
for t >tQ 

||Au(t)||^2(s.2) < (^ + (^ + (l + , (5.6) 

for some positive non-dimensional constant c. 

Remark 5.1. One can also follow the argument in [10] using the analyticity of the 
solution u{t), in time, and the Cauchy integral formula to derive an estimate on 
||^||^2(Q) • The estimate on ||Am||^ 2 ^qj will then follow by a similar argument as 
presented above and bounded by a similar estimate as in (5.6). 
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